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i!B ■ ABSTRACT 

Ph. We present the general formalism for studying CP-violating phenomena in the production, 
• mixing and decay of a coupled system of CP-violating neutral Higgs bosons at high-energy 
colliders. Considering the Minimal Supersymmetric Standard Model (MSSM) Higgs sector 
^ ' in which CP violation is radiatively induced by phases in the soft supersymmetry-breaking 
5y I third-generation trilinear squark couplings and gaugino masses, we apply our formalism 
to neutral Higgs production via bb, gg and W~^W~ collisions at the LHC. We discuss CP 
asymmetries in the longitudinal and transverse polarizations of r+r" pairs. The signatures 
of CP violation are more prominent in the production via gg and W^W~ than via bb, and 
are resonantly enhanced when two (or all three) neutral Higgs bosons are nearly degenerate 
with mass differences comparable to their decay widths. Such scenarios occur naturally in 
the MSSM for values of tan/3 > 5 (30) and large (small) charged Higgs-boson masses. 
We analyze representative examples with large mixing between the three neutral Higgs 
bosons weighing about 120 GeV, that may exhibit observable CP asymmetries even as 
large as 80%. 
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1 Introduction 



If super symmetry (SUSY) turns out to be realized at low energies < 1 TeV [1], the following 
interesting questions will then arise: does SUSY make observable contributions to the 
violation of either flavour or CP? Even in the minimal supersymmetric extension of the 
Standard Model (MSSM), the soft SUSY-breaking sector may include about a hundred 
parameters that violate these symmetries. However, if one imposes flavour conservation on 
the soft SUSY-breaking parameters mo,mi/2 and A, and assumes that they are universal, 
then only two physical CP-violating phases remain: one in the gaugino masses mi/2 and 
one in the trilinear couplings A. 

These CP-violating phases may in principle be measured directly in the production 
cross sections and decay widths of sparticles at high-energy coUiders [2,3], or indirectly via 
their radiative effect on the Higgs sector [4] ^. The Higgs sector of the MSSM is affected at 
the one-loop level by the trilinear phase [4, 10-17] and at the two-loop level by the gaugino 
mass phase [12-15]. This loop- induced CP violation mixes the CP-even Higgses h, H with 
the CP-odd Higgs boson A. Many studies have been made of the masses and couphngs of the 
resulting mixed-CP Higgs bosons -f^i,2,3, and some of their phenomenological consequences 
for searches at LEP and future coUiders have also been considered [18-26]. 

More complete studies of CP-violating Higgs bosons will require a careful treatment 
of the resonant mixing of multiple Higgs bosons that couple to the same initial and final 
states. In general, one could expect that the CP-violating mixing of the heavier MSSM 
Higgs bosons H, A may be more important than their mixings with the lightest CP-even 
Higgs boson h. However, non-negligible mixing among all three neutral Higgs states is 
also possible in a general CP- violating MSSM. Such a scenario naturally emerges from a 
parameter space where tan /3 is large, i.e. tan /3 > 30, and the charged Higgs bosons are 
relatively light with Mh± ^ 160 GeV. 

In this paper, we develop the general formalism for describing the dynamics that 
governs the production, mixing and decay of a coupled system of CP- violating neutral Higgs 
bosons. Our formalism makes use of the field-theoretic resummation approach developed 
in [27] to treat unstable particle-mixing transitions. Within the context of gauge theories, 
it is important that resummation approaches to unstable particles consistently maintain 
crucial field-theoretic properties, such as gauge invariance, analyticity and unitarity [28]. 
It has been shown in [28, 29] that all these properties are preserved within the framework 
of the Pinch Technique (PT) [30]. Here, using the PT, we compute the gauge- mediated 

^Additional indirect constraints on the soft SUSY-breaking phases and the MSSM mass spectrum may be 
obtained from experimental limits on electric dipole moments (EDMs) [5-7] and B-meson observables [8,9]. 
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diagonal as well as off-diagonal absorptive parts in the resummed Higgs-boson propagator 
matrix. Finally, an essential ingredient of our formalism is the inclusion of the CP-violating 
loop corrections in the production and decay vertices of the Higgs bosons. 

We illustrate our general formalism for the coupled-channel -f^i,2,3 mixing by explicit 
treatments of the production processes gg, bb ^ and W^W~ -^^1,2,3 t^t~ . These 
are the most important production mechanisms for neutral Higgs bosons at the LHC, 
while the decay channel that seems the most promising for studies of CP violation is that 
into r+r~ pairs. To quantify the genuine signatures of CP violation, we calculate CP 
asymmetries that are defined in terms of longitudinal and transverse polarizations of the 
leptons. When tan/5 is large and/or the charged Higgs boson mass is large, so that 
two or more Higgs bosons are nearly degenerate, even small CP- violating phases could 
induce sizeable CP-violating mixing. However, as we demonstrate, there are systematic 
cancellations due to CPT-preserving rescattering effects in the process bb -^^1,2,3 ^ r+r^ 
that suppress the CP-violating signatures in this case. There are no such cancellations in 
gg and W^W~ —>■ -f^i,2,3 r^r^, which could have much larger CP asymmetries at the 
LHC. We analyze representative examples with large three-way mixing to show that these 
CP asymmetries might well exceed the 10% level and could even reach values up to 80%. 

Although our predictions are obtained in the MSSM with explicit CP violation, it is 
important to stress that large CP- violating effects could also occur in a general CP- violating 
2-Higgs-doublet model with a similar Higgs-boson mass spectrum. Our presentation is 
organized in such a way that the formalism may easily be extended to Higgs production at 
other colliders. For instance, our formalism can be applied to 77 colliders [25,32], which are 
analogous to gg collisions at the LHC, to fi^fi~ colliders [26], which have formal similarities 
with bb collisions at the LHC, and to PFiy-fusion and Higgsstrahlung processes at e'^e' 
linear colliders [33]. 

Section |21 presents the general formalism for the coupled-channel analysis of Higgs 
bosons, including explicit formulae for the absorptive parts of the Higgs-boson propagator 
matrix and the vertex corrections. In Sectional we apply the results of our formalism to the 
production channels gg, bb, W~^W~ —>■ t^t^ at the LHC. In Section|3]we present numerical 
estimates of particular CP-violating MSSM scenarios that exhibit large CP asymmetries. 
Our numerical estimates are based on the Fortran code CPsuperH [34]. Finally, Sectional 
contains our conclusions and discusses the prospects for pursuing studies of Higgs-sector 
CP violation at future colliders beyond the LHC. 

^ We note that the hb fusion process may become the leading production channel at large tan (3 at the 
LHC, as has recently been shown in [31]. 
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2 Formalism for Coupled-Channel Analyses of Higgs- 
Sector CP Violation 



We consider situations where two or more MSSM Higgs bosons contribute simultaneously to 
the production of some fermion-antifermion pair whose polarization states can be measured. 
We treat explicitly the example of -f^i,2,3 t~^t~ , but the formalism could easily be adapted 
to other cases such as tt, xtxj or x^iX^- There have been extensive discussions of the 
masses and couplings of MSSM Higgs bosons mixed by loop-induced CP violation [10-13]. 
To account properly for the constraints that CPT invariance and unitarity imposes on the 
cross sections [27], we must consider the full off-shell propagator matrix for mixed MSSM 
Higgs bosons, including off-diagonal absorptive parts. 

The absorptive part of the Higgs-boson propagator matrix receives contributions from 
loops of fermions, vector bosons, associated pairs of Higgs and vector bosons, Higgs-boson 
pairs, and sfermions: 

$>mni,(s) = QmU{l (s) + QmUY/'{s) + QmU^^{s) + QmUfj"{s) + Qmli{/{s) . (2.1) 

The contributions of the exchanges of the bottom and top quarks, r leptons, neutralinos 
x'i and charginos xt summed in '^mn{j (s). The latter may conveniently be cast into 
the form 

^mUilis) = £ Yl Kfis)g}Aff,Nl,[{l-Kf-Kf>)igyfgf;j,f + g^j,fg^^^^^^ 

f,f'=b,t,T,x<',X~ 

-2y/KjKj7{gHj>f9H*ff - 9Hj'f9H*f-'f)] '^^^^(l' i^f^ i^f) ^{s- {rrif + m^f) , (2.2) 

where K,,{s) ^ 1+5.67^, Aff, = Sfj, (/, /' = 6, t, r), ^ (/, /' = x?,2,3,4), or 1 (/, /' = 
Xi,2)j ^{x-iy-i^) = x^ + y'^ + z'^ — 2{xy + yz + zx) with = rri^/s. Here and subsequently, 
we follow the convention of CPsuperH [34] for the couphngs of the Higgs bosons to fermions, 
vector bosons, Higgs bosons, and sfermions. For the calculation of the bottom- and top- 
quark contributions, the running quark masses at the scale ^/s have been used in the 
couplings gb^t — g'nib,t{\fs) /'^Mw- Specifically, we use m;,(mf°'*') = 3 GeV, where mf°'^ = 
175 GeV. 

The vector-boson loop contributions are 

9'^9h,vv9„,vv^V 



+2M-UMI + Ml.) + MlMl^ Q{s - AM^), (2.3) 
where /3y = (1 - 4^^)^/^ and 5w ^2, 5z ^ 1. 
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Correspondingly, the exchanges of Higgs and vector boson pairs give 

D47riW^ 2,3 

+(M| - MlJiMl + Ml^) + MlMl] Q{s- {Mz + Mh,] 



J' 



+ (M^ - M^±)(M4 + M4.) + M^M^.J e[s- {Mw + Mh±Y) . (2.4) 

In deriving (j2.3|) and (|2.4|) . we apply the PT to the MSSM Higgs sector following a procedure 
very analogous to the one given in [35] for the SM Higgs sector. As a consequence, the 
PT self-energies 9mH]^^(s) and 53mHj^^(s) depend linearly on s at high energies. This 
differs crucially from the bad high-energy dependence oc that one usually encounters 
when the Higgs-boson self-energies are calculated in the unitary gauge. In fact, if the 
Higgs-boson self-energies are embedded in a truly gauge-independent quantity such as the 
S-matrix element of a 2 — 2 process, the badly high-energy-behaved s^-dependent terms 
cancel against corresponding terms present in the vertices and boxes order by order in 
perturbation theory. In this context, PT provides a self-consistent approach to extract 
those s^-dependent terms from boxes and vertices, thus giving rise to effective Higgs self- 
energies that are independent of the gauge-fixing parameter and s^. More details on the 
PT may be found in [28-30,35]. 

Finally, the contributions of the MSSM Higgs bosons and sfermions are 

(2.5) 

-L"-"' f=b,t,T k,l=l,2 

(2.6) 

Note that the symmetry factor Sij-ki has to be calculated appropriately. When i = j = 1 
and k = I = 2, for example, the symmetry factor for the squared self-coupling S'^ ^ ^ is 

'S'll;22 = 4. 

When considering any specific production process and decay channel, the Higgs-boson 
propagator matrix must be combined with the appropriate vertices, that themselves receive 
CP-violating loop corrections. Since the main decay channel we consider for the LHC is 
r+r^, and since many of the interesting Higgs production and other decay mechanisms 
also involve fermions such as bb, we also summarize relevant aspects of the loop-induced 
corrections to the -^1,2,3// vertices. 
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The exchanges of gluinos and charginos give finite loop-induced corrections to the 
Hi 2 366 Yukawa coupling with the structure 



hh 



V2mb 1 

vcosP 1 + {5hb/hb) + (Ahb/hb) tan/? ' 



(2.7) 



The tan /3-enhanced threshold correction (Ahb/hb) has terms proportional to the strong 
coupling as and the top-quark Yukawa couphng See Eqs. (2.4) and (2.5) in [20] for the 
analytic forms of (Shb/hb) and (Ahb/hb), respectively. In addition, there are contributions 
to (Ahb/hb) coming from the exchanges of binos and winos which are proportional to the 
electromagnetic fine-structure constant ctem [36]. Taking CP violation into account, these 
additional contributions read 



(Ahb/hb) 



aemlJ'*M* 



\UL\'l(7nl\M,\M^i\') + \UU'l(ml\M,\M^\ 
1 



+ ^\UL\'I(mi^,\M,\\\l^\') + - \Ul\'I(ml,\M,\\\i^\') 



ll(ml,ml,\Mi\'' 



l\ui\'i(mim^\^\^) 



+ -\Ul,\'l(ml,\Mi\',\i,\') + |C/^i|2/(m?,|Mi|M//|2) 



+\K2\'nK:m\\\pi\') 



where 



I(a, b, c) 



ab ln(a/6) -|- be ln(6/c) -|- ac ln(c/a) 



(2.8) 



(2.9) 



(a — b) (b — c) (a — c) 

We follow the convention of CPsuperH [34] for the mixing matrices of the stops sbottoms 
and staus U^. 

There are formulae analogous to those above for the loop corrections to the Hi^2,3ti 
vertices, which would be relevant for CP- violation measurements in e~e~^ uutt [33], for 
example. 

Analogous exchanges of binos and winos give finite loop-induced corrections to the 
Hi 2 3T~^T~ coupling, which have a similar structure: 



where 



(^hr/hr) 



aem/^*M2* 



4:77 s'^ 



w 



V COS (3 1 + (Ahr/ hr) tan (5 ' 



I(ml,\M2\'M') + -\Uli\'I(mi^,\M2\'M') 



(2.10) 
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\M2\\ 













+ 



4vr 



/K,m?,|Mi|2) 



-|t/^J^JK^,|Mi|M/x 



[/^2rA<,|Mi|M/x| 



(2.11) 



The threshold corrections modify the couphngs of the neutral Higgs bosons to the scalar 
and pseudoscalar fermion bilinears as follows [20]: 



1 



o 



P 



cos P 
+ Oai S^m 
O. 



1* 



cos P 



where f = h and r and 



Kb 



1 + tan P J 
/ Kf (tan2/5+ 1)' 
\ 1 + tan P 
/ Kf tan/9 \ 
\1 + K/tan/3 J 
/ tan /9 — \ 
yl + Kf tan P J ' 

(A/ib/Zib) 



cos P 



3?e 



1 + Kf tan P ^ 



O 



cos P 



1 + tan P ^ 



Kr = {Ahr/hj 



(2.12) 



(2.13) 



1 + {6h/h) ' 

There are similar formulae for the iJi.2,3/i^/i~ vertices that would be relevant for ^^fi^ 
colliders. The analogous corrections to the Hi 2^3e~^e~ vertices may be neglected. 

Additional contributions to Higgs-boson vertices may arise from absorptive effects 
due to the opening of third-generation sfermion pair production channels. However, if the 
-f^i.2,3-boson masses are well below the kinematic threshold of such production channels, 
these absorptive effects are small and can be neglected. Finally, we remind the reader 
that detailed analytic expressions for the effective Higgs-boson couplings to the photon, 
the gluon, the bosons and SUSY particles are given in [34]. 



3 Tau Pair Production at the LHC 

To further elucidate the formalism presented in Section |2l we now discuss in more detail 
the production, mixing and decay of Higgs bosons into polarized r+r^ pairs at the LHC. 
We will study individually the three most significant production channels for Higgs bosons 
in the MSSM at the LHC: (i) bb fusion, (ii) gg fusion and (iii) W^W~ fusion. 
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Figure 1: Mechanisms contributing to the process bb H ^ t^t , including off-diagonal 
absorptive parts in the Higgs-boson propagator matrix. 

3.1 bb Fusion 

At large tan (3, an important mectianisni for producing neutral Higgs bosons at tlie LHC 
is bb fusion [31,37-39]. Figure illustrates how the matrix element A^^* for bb —>■ r^r^ 
receives contributions from the s-channel exchanges of the neutral Higgs bosons. The loop- 
corrected propagator matrix and vertices calculated in the previous section are indicated 
by shaded circles. The matrix element can be written as 

"^^^^W^ i,j=l,2,3 a,l3=± I 

X(^i,r+r- + iP9H,r+r-) u{pi, Ct) Ppv{p2, a)^ , (3.1) 

where = (1 + q;75)/2 and the running bottom-quark mass at the scale of is used. 
We denote the helicities of r~ and r+ by a and a and those of the b and 6 by A and 
A, respectively, with a, A = + and — standing for right- and left-handed particles. The 
four-momenta of the r~ and r"*" are pi and p2, respectively, those of the b and b are ki and 
k2, respectively, and s is the centre-of-mass energy squared of the bb pair that fuses into a 
Higgs boson: s = {ki + ^2)^ = {pi +^2)^- 

An important element of our formalism is the consideration of the 'full' 3x3 Higgs- 
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Figure 2: The t^t~ production plane with definitions of the scattering angle 0. The 
transverse polarization vectors Pt and Pt have azimuthal angles a and a, respectively, 
with respect to the event plane. 

boson propagator matrix D{s) in ()3.1|) ^. This is given by 

m = 



-1 



f s~ Ml^+i^mRii{s) i^mRi2{s) i^YnIii^{s) \ 

§ iS5mn2i(s) s- Ml^ + i'^mR22{s) i'^mli2^{s) 

\ i$jmn3i(s) iSjmn32(s) s-M^3 +iSjmn33(s) / 

(3.2) 

where the inversion of the 3-by-3 matrix is carried out analytically. In ()3.2|) . the absorptive 
parts of the Higgs self-energies 55mnjj(s) are given in Section |2l and Mh^ ^ .^ are the one- 
loop Higgs-boson pole masses, where higher-order absorptive effects on Mh-^ 2 3 have been 
ignored [15]. In the same context, the off-shell dispersive parts of the Higgs-boson self- 
energies in the Higgs-boson propagator matrix D{s) have also been neglected, since these 
are formally higher-order effects and very small in the relevant Higgs-boson resonant region. 
Finally, we include in (|3.1|) the finite loop-induced corrections to the couplings of Higgs 
bosons to h quarks, fi'^^j,, and r leptons, g^^^+^., due to the exchanges of gauginos and 
Higgsinos, as has been discussed in Section |21 

In the centre-of-mass coordinate system for the bb pair, the helicity amplitudes are 
given by 

M'\aa;XX) = (a; A),<5.Aa , (3-3) 



where 

i,j=l,2,3 

•^Strictly speaking, the complete propagator matrix D{s) is a 4 x 4-dimensional matrix spanned by the 
basis {Hi, H2, H^jG^) [27]. However, to a good approximation, we may neglect the small off-resonant 
self-energy transitions of the Higgs bosons Hi 2.3 to the neutral would-be Goldstone boson G°. 
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with i3f = y 1 — Am'j/s. Note that the cross sections for general (longitudinal or transverse) 
polarizations can be computed from the helicity amplitudes by a suitable rotation [40] 
from the helicity basis to a general spin basis. 

The r-polarization weighted squared matrix elements are given by 

1 



M 



bb 



12 



12 



E Tr M''fM''^p 



(3.5) 



A=± Wcto-' / A=± 

where p and p are 2x2 polarization density matrices for the r~ and r"*", respectively: 



P 



1 + Pl Pre- 



1 



1 + 



1-Pl 

Evaluating the trace in ()3.5|) yields 



P 



2 I _p. 



re 



1-Pl 



(3.6) 



7W 



66 



g mbm-j- \ 



12 I 



C?(l + PiPi) + C2^(Pz. + Pi) 



+PtPt C\ cos(a - a) + sin(a - a) 



(3.7) 



The T'^T production plane is depicted schematically in Fig. |2l where the transverse po- 
larization angles a and a are also defined. 

The coefficients C^(n = 1 - 4) in (jSHI) are defined in terms of the helicity amplitudes 

by 



^ A=± 



cl ^ ^E (K+;^)^l'-l(-;^)^ 



A=± 



Cl 



i E ((+; A),(-; A)^) , ^ ^ E ((+; A),(-; A)^) . (3.8) 



A=± ^ A=± 

Under CP and CPT ^ transformations, the helicity amplitudes transform as follows: 



(or; A), S -\), , {a; X)t "A^ +{-cr; -A),* . (3.9) 

Hence, the CP and CPT parities of the coefficients defined in ()3.8|) are given by 



(3.10) 



where the first and second symbols in the square brackets are the CP and CPT parities, 
respectively. Consequently, the coefficients C\ and C\ signify genuine phenomena of CP 



"^We define T as the naive T-reversal transformation, under which the spins and 3-momenta of the 
asymptotic states reverse sign, without interchanging initial to final states. In addition, under the operation 
of T, the matrix element gets complex conjugated. 
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violation, whereas and C| are CP-conserving. Here, we should remark that a non-zero 
value for the CPT-odd coefficient can only be induced by non-vanishing absorptive 
effects. In our case, such effects mainly originate from the absorptive parts of the Higgs- 
boson self-energies. 

Finally, for our phenomenological discussion in Section 01 we define the parton-level 
cross sections 

where the intermediate state H collectively denotes all the Hi Hj resonant transitions 
with i,j = 1, 2, 3. 



3.2 gg Fusion 

The matrix element Ai^^ for the process gg H t^t^ , depicted in Fig. El can be 
written as 

SnvMw 

3 

^ II GnXh, ei; ^2, £2) Ai(s) {gH.r+r- + ^(^9h,t+t-) (T)PaV{p2, ^) ■ (3.12) 

i,j=l,2,3 a=± 

In the above, a and h are indices of the SU(3) generators in the adjoint representation and 
ki^2 and ei^2 are the four-momenta and wave vectors of the two gluons, respectively. Again, 
we denote the helicities of r~ and r"*" by a and a with a = + and — standing for right- 
and left-handed particles. The four- momenta of t~ and r+ are pi and p2, respectively, and 
s = (/ci + ^2)^ = (yPi + P2Y- The Higgs-boson propagator matrix D{s) was given in ()3.2j) 
and the loop-induced couplings of the Higgs bosons Hi to two gluons are given by 

GHXki,(^uk2,e2) = i Sf {y/I) (^ei ■ e2 - ki ■ e2k2 ■ ei^ ~ i Pf{\/§) e^^p^e'^e'^k'^k^ , 

(3.13) 

with £0123 = 1- For the loop functions Sf and Pf, we follow the definitions of [34]. 

In the two-gluon centre-of-mass coordinate system with ki along the positive z di- 
rection and k2 along the negative z direction, the wave vectors of two photons are given 
by 

e^'(Ai) = -^(0,-Ai,-«,0), e^(A2) = (o, -A2, 2, o) , (3.14) 
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where A = +1 and —1 denote the right and left gluon hehcities, respectively. The helicity 
amplitudes are given by 

M^\aa;\,\,) = ^^f^l^^{a; \,) ,6^,5,,,, , (3.15) 

oT[VlVl\Y 

where the amplitude (a; X)g is defined as 

(a;A),^ E [S!{VI)+^XPf{VI)]D,,{s){cr(3^g|,^,,,--^g^^,,,-). (3.16) 

i,j=l,2,3 

We note that the amplitude ()3.15p has the same structure as the amplitude ()3.3|) for 
bb T^T~, except for the overall constant. We obtain from the helicity amplitudes the 
polarization-weighted squared matrix elements given by 

32 Y 8ttvM\y J [ 

+PtPt [CI cos(a -a) + Ci sm{a - a)] | , (3.17) 

where the coefficients are obtained by replacing (cr; \)b ^ (cr; A)^ and interpreting A as 
the gluon helicity in (j3.8p . 

Under the CP and CPT transformations, the helicity amplitudes transform as follows: 

(cr; A), S -{-a; -A), , (a; A), W -{-a; -X); , (3.18) 

where the CP and CPT parities of the coefficients Cf are the same as those of the 
Finally, we define the parton-level cross sections as: 

a.igg^H^r^r-) ^ ^ f ^^f!^ V . (3.19) 
' 512ns V SnvMw J 

12 



Note that the CP- and CPT-odd cross section (T2 receives contributions from the absorptive 
parts of the Hi^2,3gg vertices and Higgs-boson self-energies as well. 

3.3 W+W- Fusion 

The last important mechanism for the production of the MSSM neutral Higgs bosons at the 
LHC is W^W- fusion [41-44]. The matrix element M^^ for this process, W-{ki)W+{k2) - 
H T~ {pi)T^ {p2) with s = {ki + k2y, is given by 

^WW ^ g^^^^ . {9H.r+T- + i<^9HjT+T-) U{P1, (T)PaV{p2: a) , 

i,j=l a=± 

(3.20) 

where ei and 62 are the polarization vectors of two vector bosons and Qh^vv denotes the cou- 
pling of the Higgs boson Hi with a pair of gauge bosons, as defined through the interaction 
Lagrangian 

Chvv = gMw iw^W-^' + ^Z^Z'^\ QH^vvHi . (3.21) 

In the W^W'' centre-of-mass coordinate system with ki along the positive z direction and 
k2 along the negative z direction, the polarization vectors of two vector bosons are given 
by 

e5'(Ai = ±l) = i^(0,Tl,-i,0), ei'(Ai = 0) = ^ (|ki|, 0, 0, A;?) , 
v2 Jk\ 



e^(A2 = ±l) = ^(0,:Fl,i,0), e^(A2 = 0) = (|k2|, 0, 0, , (3.22) 

V2 ^kl 

where the polarization vectors are normahzed by ei(A) •e*(A') = — ^aas and A = ±1 and A — 
denote the transverse (right and left helicities) and longitudinal polarizations, respectively. 
In this frame, the helicity amplitude is given by 

2 

M^^((ja;AiA2) = ^-^^a-M)w?>c-.K^,. (3.23) 

2vs 

where the amplitude (cr; A)vk is defined by 

(a; \)w = I] '^(A) gH.vv Ai(s) [p^r 9%r+T- - Whjt+t-) > (3-24) 

i,j=l,2,3 

with 



uj{±) = 1 and a;(0) = -ki ■ k2/Jklkl . (3.25) 
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The factor ci;(0) becomes 1 — s/2M^ for on-shell vector bosons and dominates the amphtude 
for s > M^. 

One can then obtain the following averaged amplitude squared: 



\M 



+PtPt cos(a -a) + Cf sin(a - a) 



(3.26) 



where the coefficients C„ can be obtained by replacing (a; A);, {a; X)w and summing 
over A = ±,0 in ()3.8j) . The CP and CPT parities of the coefficients are the same as 
those of or C^, and the parton- level cross sections are defined similarly as 



/3r 



144 TTs V 2v/I 



C. 



w 



(3.27) 



In kinematic situations where the longitudinal W^W£ contributions can be neglected, the 
average factor 1/144 should be replaced by 1/64. Finally, we note that it is straightforward 
to calculate ZZ-fusion processes in a similar fashion, although their cross sections are 
smaller approximately by a factor of 4 than W^W~ collisions at the hadron level. 



4 Numerical Examples 

We now present some numerical examples of CP-violating Higgs signatures in t^t~ produc- 
tion at the LHC. As already mentioned, these signatures may be enhanced at large tan/3, 
and three-way mixing is potentially important for small charged Higgs-boson masses. Since 
the prospects for observing -^1,2, 3 t^t~ at the LHC are best for light Higgs bosons, we 
present in this section some numerical analyses in a specific scenario in which all the three 
Higgs states mix significantly. 

Explicitly, we take the following parameter set: 

tan/3 = 50, MP±" = 155 GeV, 

= Mfj^ = = = = MsusY = 0.5 TeV, 
= 0.5 TeV, IAatI = 1 TeV, IM2I = |Mi| = 0.3 TeV, IM3I = 1 TeV, 
$^ = 0°, $^ = $^^=$^^ = $^^=90°, $i = $2 = 0°, (4.1) 

and we consider two values for the phase of the gluino mass parameter M3: $3 = —90° , —10°. 
For $3 = —10°, CPsuperH yields for the masses and widths of the neutral Higgs bosons: 
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Mh, = 120.2 GeV, Mh, = 121.4 GeV, Mh, = 124.5 GeV, 

Th, = 1.19 GeV, Th, = 3.42 GeV, Th, = 3.20 GeV, (4.2) 

and for $3 = -90°: 

Mh, = 118.4 GeV, Mh^ = 119.0 GeV, M^g = 122.5 GeV, 

Th, = 3.91 GeV, = 6.02 GeV, Th^ = 6.34 GeV, (4.3) 

respectively. 

In Figs. 0] and we show the parton- level cross sections (Tj(66 —>■ H —>■ r+r"), 
a^igg ^ H ^ r+r") and ai{WW H ^ r+r") defined in flTTTIl . STI^ and (ITTfll . 
respectively, as functions of the r+r~ invariant mass V^. The solid lines are for $3 = —90° 
and the dashed (red) ones for $3 = —10°. We recall that non- vanishing of 0^2 and 
are direct signals of CP violation in longitudinally and transversally polarized r+r~ pairs, 
respectively. 

The parton- level cross sections (Ji(WW — * r+r~) have been computed by neglecting 
the contribution of the longitudinally-polarized W^^ i.e., setting ci;(0) = 0. For the MSSM 
scenario defined in (jHH), this is a plausible approximation for Higgs-boson masses below 
the WW threshold. Possible uncertainties that such a treatment may introduce largely 
cancel when we consider ratios of the cross sections aiiWW t^t^), such as the CP 
asymmetries to be defined later in this section. 

In Fig. 0J we observe that the cross section (3"2, which quantifies CP violation in the 
production of longitudinally polarized r-lepton pairs at the parton level, is comparable to 
the spin-averaged cross section ai in WW and gg collisions. This implies that CP violation 
can be very large in these channels. Instead, in bb fusion, the ratio (T2/(Ji is always less 
than 1%, so CP-violating effects in the production of longitudinally polarized r leptons are 
unobservably small in this case. 

The smallness of (J2 in bb fusion is a result of an intriguing interplay between unitarity 
and CPT invariance [27]. In detail, the CP- violating cross section o"2 may be calculated by 

1 



a2{bb ^ r+r-) ^ 



a(bb ^H^ r+r^) - a{bb ^ H r^r^ 



(4.4) 



where a denotes the usual subprocess cross section. For the scenario under study, unitarity 
cuts of bb pairs dominate the absorptive part of the Higgs-boson self-energies. Employing 
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this fact and the optical theorem, we obtain 



\=L,R ^ 

- ^mTir+TE ^H^ r+r^) ) + 0\a',B{H,,2,3 ^ r+r")] , (4.5) 



where Cps is a phase-space correction factor and T{t^j^Ti^j^ H ^ T~L,R^L,ii) denote 
the usual matrix elements. In ()4.5|1 . a'2 is the CP- violating cross-section (3"2 calculated by 
omitting the off-diagonal absorptive parts in the Higgs-boson propagator matrix D{s). The 
size of (3"2 is smaller at least by a factor 10 than the spin-averaged cross-section ai. On the 
other hand, CPT invariance imposes the constraint 

T(r+r^ ^H^ r+r^) = T(r+r^- ^H^ rtrl) . (4.6) 

With the aid of ()4.6|) . it is not difficult to see using ()4.5p that the CP- violating cross section 
(72(66 H T^T~) vanishes up to CP- violating terms suppressed by extra factors of order 
5(i/i,2,3 ^ r+T-). 

Our numerical estimates presented in Fig. El show that the CP- violating transverse- 
polarization cross section may be quite sizeable for all production channels. However, 
generically exhibits an alternating sign for hb and gg collisions, and CP violation becomes 
very small when we integrate over the whole Higgs-boson resonance region. Moreover, 
the transverse polarizations will be difficult to measure at the LHC because of the 
experimental conditions, notably the large boosts of the r^. On the other hand, analogous 
asymmetries might be observable in it production and/or in t^t~ production at a 
collider. 

The Higgs production channels via bb and gg fusion processes can be separated from 
the W^W~ fusion channel by applying a number of kinematic cuts [45] including the 
imposition of a veto on any hadronic activity between jets [46,47]. Therefore, we treat 
the contributions from hb and gg collisions to the physical Higgs-exchange process pp —>■ 
H T~^T~X separately from those coming from WW fusion. More explicitly, the physical 
r+r^ cross section can be computed by integrating the parton- level cross sections with the 
distribution of b quarks, gluons and VT-bosons in the proton, 

r^(^pp{bb,gg)^T+T-x) = Aa,{bb ^ H ^ r+r-) r ^ 

+ AKa,{gg^H^ r+r-) r , (4.7) 

r 7 rWW 

r'^ippiW+W-) ^ T+T-X) = AaiiW+W- ^ H ^r+r-) r—^ ,(4.8) 
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where r is the Drell-Yan variable t = s/s and s is the invariant squared centre-of-mass 
energy of the LHC. In ()4.7|1 . we use the value K = 1 + (vr^ + 11/2), ignoring the small 
difference between the fC— factors for CP-even and CP-odd Higgs states. The effective 
luminosities for bb and gg collisions, and £^^, may be determined by 



b{x,Q)b(^^,Qj + {b^b) 



r — ; — = ax 
dr Jt I X 

= j^dx-g{x,Q)g(^-,Qy (4.9) 

where b{x,Q), b{x,Q) and g{x,Q) are the b, b and gluon distribution functions in the 
proton and Q is the factorization scale. In our numerical analysis, we use the leading-order 
CTEQ6L [48] parton distribution functions for b{x,Q) and b{x,Q), and the CTEQ6M 
parton distribution function for g{x,Q). We choose the factorization scale Q = -\/^/4 for 
the 6-quark fusion process as suggested and confirmed in [38]. 

Correspondingly, in ()4.8p . the effective luminosities for the transverse and longitudinal 
ly-bosons, denoted as W^^^, can be computed in terms of effective densities ± (x, Q) in 
the colliding protons, which are in turn calculated in terms of the quark parton distribution 
functions q{x, Q) in the proton:^ 



(l£WpWp ^1 

r 



dr 



F'w-M) = E f'-<liy,Q)Fl4^/y,Q), (4.10) 



g=u,d,c,s 



T,L 



where the transverse (P = T) and longitudinal (P = L) effective densities in the 

quark q are given by [42,43] 

J-^n I ^\ Q^em , i Q \ 1 "I" (1 



X 



K-M) = ^l—^. (4.11) 

4:71 X 



Note that the summation over quark flavours q in the expression for (x, Q) includes 

''T,L 

q = u,d, c, s. Moreover, we take Q = in our numerical estimates. 

To analyze the signatures of CP violation in the production of longitudinally polarized 
r-leptons, we first define the physical observables 

(^RR = a{pp ^ H ^ r^rgX), (JLL = (t{pp ^ H ^ rlr^X). (4.12) 



^Here we consider identical polarizations for the bosons in the W'^W fusion process. 



17 



Evidently, the total cross section for Higgs production and decay into r+r pairs is given 
in terms of a^R and ctll by 

CTtotipp ^ H r+r"X) = aRR + gll ■ (4.13) 

Although the initial state pp is not symmetric under CP, it can, however, be shown that, 
up to negligible higher-order CP-violating electroweak effects, the effective luminosities for 
gg^ bb and W^W~ densities will be practically the same for pp and pp collisions. Therefore, 
the difference of cross sections 

A(TCP = (JRR - ctll (4.14) 

is a measure of genuine CP violation at the LHC. In analogy with ()4.7|1 and ()4.8j) . the 
CP-violating cross section Aacp can be computed by 



bb 



dr ^ ' dr 

dC^^ 

+ AKa2{gg^H^ r+r") r , (4.15) 

dArr ^rWW 

T=^(pp{W+W') T+T-X) = 4a2(W+W- ^ H ^ r+r-) T—, . (4.16) 

dr ^ ' dr 

To gauge the sizes of the signatures of CP violation at the LHC, we define the following 

two CP asymmetries: 

^ dr A - ^^cp 



acp(r) ^ fj , ^cp = , (4.17) 

r 



(Mot (Ttot 

dr 



pertinent to the hadron- level processes (66, (7(7, VTW^) ^ H ^ t^t X. 

We plot in Fig. IHl the differential cross sections r ^^^^ and as functions of 

V^. The upper two frames are for the process 66 — > r+r~, the frames in the middle for 
gg T^T~, and the lower ones for W^W~ t^t^ . We observe that the the main 
production mechanism is 66 fusion, which gives a cross section about five times larger 
than that due to gluon fusion for the scenario under consideration. However, as has been 
mentioned above, the iy"'"iy~-fusion production mechanism, albeit much smaller, can be 
experimentally distinguished from that due to 66 and gg collisions. Therefore, in Fig. Owe 
display the CP asymmetry acp defined in ()4.17|1 separately for the 66 -|- gg and W^W~ 
subprocesses. We note that the large CP asymmetry in gg subprocess is diluted by the 
dominant cross section via 66 fusion®. 



^Specifically, the total CP asymmetry in the gg subprocess is A^p = —8.4 (—6.2) %, for $3 = 
-90° (—10°). However, after the inclusion of bb collisions, the combined CP asymmetry J^^^^^ reduces to 
-1.4 (-1.0)% . 
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It is important to emphasize here that the CP-violating observable r '^^^^^ for the 
IWT-fusion process does not change sign as the r'*'r~-system energy varies over the 
entire Higgs-boson resonant region. Such a kinematic behaviour is ensured by the presence 
of the off-diagonal absorptive parts of the Higgs-boson self-energies. Instead, if these off- 
diagonal absorptive parts are neglected, we find the erroneous result that the CP-violating 
observable flips sign in the resonant region, thereby leading to unobservably small CP 
asymmetries when averaged over the energy v^- 

Although CP violation in the WW and gg production channels may be sizeable, it 
is difficult to measure the differential CP asymmetry acp at the LHC because of the low 
energy resolution of the reconstructed r+r~ invariant mass. This last fact also limits our 
ability to reconstruct with sufficient accuracy the line shape of the decaying coupled Higgs- 
boson system at the LHC. This is unfortunate since one would miss the very interesting 
feature shown in Fig. IHl that, unlike the case of a single resonance, the locations of the 
various maxima in the resonant line shapes described by r crucially depend on the 
production and decay channels of the coupled Higgs-boson system. Therefore, the extra 
analyzing power of e^e~ and ^'^ ^~ colliders would be highly valuable for unravelling the 
existence of a strongly-mixed Higgs-boson system and studying in more detail its dynamical 
properties. 

Motivated by the large differential CP asymmetry in the W^"*"14^~-fusion process, we 
perform a numerical analysis of the total CP asymmetry for the reaction pp{WW) H 
T^T^, integrated over the Higgs resonance peaks. We present in Figs. |Hl IHl and fTUl the 
predicted values for the cross-section <Jtot{pp{WW) — > /J — > r+r^X) and its associated 
total integrated CP asymmetry A^p^ defined in fl4.17|) as functions of $a = ^At = = 

for $3 = —10°, —70°, and —90°, respectively. In the upper two frames of the figures, 
we display the dependence of the Higgs-boson masses and their decay widths on the CP- 
violating phase $A, where the solid, dashed and dotted lines refer to the Hi, H2 and H3 
bosons, respectively. In our numerical analysis, we fix the remaining parameters of the 
MSSM as in fj4.1|) . Unlike in Figs. |H1 El and^J we present in Fig. ITT] numerical estimates 
by fixing the value of $a to —90°, but varying the CP- violating phase $3. For the scenario 
under study, all three Higgs bosons mix among themselves significantly, giving rise to level 
crossings as the CP-odd phases vary. These effects of level crossing lead to a non-trivial 
behaviour in Th,, which is between 1 GeV and 10 GeV and in ^J^p^. We find that the 
total cross section is between 0.1 pb and 0.7 pb and is comparable to the corresponding SM 
cross section 0.3 pb for Mh^-^ = 120 GeV [45]. We observe that the CP asymmetry A^p^ is 

"^In Fig.ini the widths of the Hi and H2 become larger than 10 GeV when $3 > 100° or $3 < -140°, 
where Mh^ decreases very rapidly and Hi decouples from the H2 — Hg, mixing system. 
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large for a wide range of CP phases and can even be as large as 80% for $3 = —70°. Even 
for small CP- violating phases, $3 = —10° and (180° — |$a|) < 20°, the CP asymmetry can 
be ~ 50%, as shown in Fig. |H1 Again, we note that possible uncertainties in the calculation 
of the cross sections largely cancel in the CP asymmetry ^cp- 

Finally, we comment briefly on the possible impact of low-energy constraints on the 
CP asymmetries, especially those arising from the non-observation of the electron and 
neutron EDMs and the absence of the Higgs-mediated 5-meson decay Bg^d fJ^fJ' at the 
Tevatron [49]. The EDM constraints may be considerably relaxed if we consider scenarios 
with the first two generations of squarks heavier than about 3 TeV, and if we allow some 
degree of cancellations [6] between the one- and higher- loop EDM contributions [50]. For 
the scenarios under study, we have estimated that the required degree of cancellation is 
always smaller than 80%, where 100% corresponds to complete cancellation. Therefore, a 
full implementation of EDM constraints will not alter the results of the present analysis 
in a significant way. On the other hand, the lack of observation of 5^ ^ — > /i/i at the 
Tevatron [49] imposes further constraints on the parameters of the CP-violating MSSM. 
However, the derived constraints are highly flavour- dependent and can be dramatically 
relaxed for certain choices of the soft SUSY-breaking mass spectrum that enable unitarity 
cancellations in the flavour space. For a detailed study, see [9]. 



5 Conclusions and Prospects 

We have presented the general formalism for analyzing CP-violating phenomena in the 
production, mixing and decay of a coupled system of multiple CP-violating neutral Higgs 
bosons. Our formalism, which is developed from [27], can be applied to models with 
an extended CP-violating Higgs sector, including the highly predictive framework of the 
MSSM with radiative Higgs-sector CP violation. An important element of the formalism is 
the consideration of the full s-dependent 3x3 Higgs-boson propagator matrix, where the 
gauge-mediated contributions to self-energies have been calculated in the framework of the 
Pinch Technique [30,35]. 

As an application of our formalism, we have studied in detail the production of CP- 
violating MSSM -ffi,2,3 bosons via bb, gg and W^W~ collisions and their subsequent decays 
into T^T~ pairs at the LHC. In addition to the Higgs self-energy effects, we have also 
given explicitly the relevant formulae in the MSSM with loop-induced CP violation in 
the production and decay vertices. We have considered speciflc MSSM scenarios that 
predict three nearly degenerate, strongly-mixed Higgs bosons with M//^ 2 3 ~ 120 GeV. 
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Such scenarios naturally occur in a general CP-violating MSSM when tan (3 is larger than 
30 and the charged Higgs boson is lighter than about 160 GeV. 

We have analyzed CP asymmetries in both longitudinally- and transversely-polarized 
T^T^ pairs. CP asymmetries that make use of the transverse polarization of the r-lepton, 
although being intrinsically very large in the CP-violating MSSM scenarios mentioned 
above, generically exhibit an alternating sign and become unobservably small after aver- 
aging over the entire Higgs-boson resonant region. Also, reconstruction of transversely 
polarized r leptons appears rather difficult at the LHC. However, such CP asymmetries 
might ideally be tested at a n'^ ijT collider, where a high energy resolution can be achieved. 

At the LHC, more promising are CP asymmetries based on the longitudinal r-lepton 
polarization. In particular, the CP asymmetry in the production channel W^W^ — > 
Hi,2,3 T^T" may well exceed the 10% level and reach values up to 80%. It is important 
to stress again that the WW production channel can be cleanly isolated from the gg and 
hh channels, mainly by vetoing any hadronic activity between jets (for details, see [45]). 
Hence, depending on the efficiency of longitudinal r-lepton polarization techniques [51], 
the production channel W^W^ -^1,2,3 t^t~ may become the 'golden' channel for 
studying signatures of Higgs-sector CP violation at the LHC. 

The formalism presented in this paper may easily be applied to other colliders as 
well, most notably to e+e^, 77 and fi^fi^ colliders. At e~^e~ linear colliders, Higgs bosons 
can copiously be produced via the Higgsstrahlung or W^W~ fusion processes. At 77 
and fi~^fi~ colliders, the polarizations of the colliding beams may also be varied, thereby 
providing additional probes of Higgs-sector CP violation. The aforementioned colliders 
can provide cleaner experimental conditions than the LHC. Consequently, even if the CP 
asymmetries discussed here prove difficult to observe at the LHC, the formalism and the 
analysis techniques developed here to investigate Higgs-sector CP violation will be directly 
applicable to such future colliders as well. 
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Figure 4: The parton-level cross sections (Ti^2{^ ^ H ^ t^t^) in pb, 0^1^2(99 ^ H ^ 
T^T^) in fb, and ai^2{W^W^ —^H^ r+r^) in pb as functions of ^/I. The solid lines are 
for the three-Higgs mixing scenario with $3 = —90° and the dashed ones with $3 = —10°. 
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Figure 5: The parton-level cross sections Cs^i^bb ^ H ^ t^t^) in pb, ^z^idd ^ H ^ 
T^T^) in fb, and (Jz^4{W~^W^ —^H^ t^t~) in pb as functions of ^/I. The solid lines are 
for the three-Higgs mixing scenario with $3 = —90° and the dashed ones with $3 = —10°. 
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Figure 6: The differential cross sections t and 

^d^cp as functions of VI. The upper 
frames are for the process bb ^ H ^ t^t^ , the middle ones for gg ^ H ^ t~^t~ and 
the lower ones for W'^W' ^ H ^ t^t^ . The solid lines are for the three-Higgs mixing 
scenario with $3 = —90° and the dashed ones with $3 = —10°. 
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Figure 7: Numerical estimates of differential CP asymmetries acp defined in as 
functions of The solid line corresponds to the three-Higgs mixing scenario with $3 = 
-90° and the dashed one to $3 = -10°. 
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Figure 8: Numerical estimates of Higgs-boson masses and decay widths, the cross-section 
atot{pp{WW) H ^ r+r^X) and its associated total CP asymmetry defined 
in as functions of = ^At = = ^Ar, for $3 = —10°. In the upper two 

frames, the solid, dashed and dotted lines refer to the Hi, H2 and H3 bosons, respectively. 
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Figure 11: Numerical values for M^^^.^ and Tu^^^, (Jx,ot{ppiW^) ^ H ^ t^t^X) and 
JiQ^ as functions of (^3, for $^ = = = ^Ar = —90°. We follow the same line 
conventions as in Fig. [3 
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